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Numerical studies on Mott transitions caused by the control of the ratio between bandwidth 
and electron-electron interaction (U) are reported. By using the recently proposed path-integral 
renormalization group(PIRG) algorithm, physical properties near the transitions in the ground 
state of two-dimensional half-filled models with the nearest and the next-nearest neighbor 
transfers (— t and t', respectively) are studied as a prototype of geometrically frustrated system. 
The nature of the bandwidth-control transitions shows sharp contrast with that of the filling- 
control transitions: First, the metal-insulator and magnetic transitions are separated each 
other and the metal-insulator (MI) transition occurs at smaller U, although the both transition 
interactions U increase with increasing t' . Both transitions do not contradict the first-order 
transitions for smaller t' /t while the MI transitions become continuous type accompanied by 
emergence of unusual metallic phase near the transition for large t'/t. A nonmagnetic insulator 
phase is stabilized between MI and AF transitions. The region of the nonmagnetic insulator 
becomes wider with increasing t'/t. The phase diagram naturally connects two qualitatively 
different limits, namely the Hartree-Fock results at small t'/t and speculations in the strong 
coupling Heisenberg limit. 
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§1. Introduction 

The Mott transition between metallic and insulator phases has been a subject of general interest 
since its proposal in the middle of the last century. The original issue was related with the failure 
of the band theory in explaining the Mott insulators such as NiO. According to the band theory, 
the highest filled band is completely filled for insulators and is partially filled for metals. The Mott 
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insulators, however, have partially filled bands and also have the insulating nature. To understand 
the insulating nature, Mott introduced the concept of electron localization caused by the mutual 
correlation effects. 1,2 ) In strongly correlated electron systems, the interplay of strong quantum 
fluctuations and correlation effects has an important role in determining the nature of materials. 
Additional complexity arises even for single-orbital systems in connection with the spin degrees 
of freedom since the electron localization itself leaves the spin entropy finite. Slater considered 3 ) 
the possibility of antiferromagnetic (AF) symmetry breaking to release this entropy, while the AF 
order itself may become the origin of the insulating gap formation. Then the metal-insulator and 
AF transitions could be coupled each other. Although the Mott transition is one of the typical 
quantum phase transitions caused by this interplay and many studies have been done, 1 ) the physical 
properties near the Mott transition have not been sufficiently understood in the presence of strong 
fluctuations. 

For the Mott transitions, there are three basic control parameters, band filling, bandwidth and 
dimensionality. 1 ) Although the quantum Monte Carlo (QMC) results in the two-dimensional Hub- 
bard model have contributed in clarifying the nature of the filling-control metal-insulator transi- 
tion, 4 ~ 6 ) the method generically suffers from the negative sign problem, if the longer-ranged transfer 
in the model becomes large. Recently we have developed an algorithm called path-integral renor- 
malization group (PIRG) method, 7 ' 10 ) which can be applied to any lattice structure without the 
sign problem. It has opened the possibility of studying correlated electron models in a much wider 
class. In this paper we apply PIRG method for studying the Mott transition which has not been 
able to be studied in other numerical methods. In particular, we study metal-insulator and mag- 
netic transitions, in two-dimensional systems, caused by the control of the bandwidth relative to 
the strength of the on-site electron-electron interaction U. The simplest Hubbard model on the 
square lattice with the nearest-neighbor transfer t satisfies the perfect nesting condition and is 
expected to be the antiferromagnetic insulator at any nonzero U/t. In real materials, however, the 
longer-ranged transfer is generically present in tight-binding representations, which drives a metal- 
insulator transition at nonzero U/t. To clarify the nature of the bandwidth-control transitions, we 
study the phase diagram in the parameter space of the bandwidth relative to U and the strength 
of the next-nearest neighbor transfer t'/t. 

Systems with the two-dimensional anisotropy are seen in many materials such as high-T c cuprate 
superconductors, organic compounds, some materials under strong electronic fields and so on. In 
more general, the bandwidth-control transition is also seen in many other systems such as V2O3 
and RMO3. 1 ) In this study, however, we do not focus our interest on any concrete materials and 
rather employ two-dimensional single-band Hubbard models as our target to study the general 
physics of quantum phase transitions in electron systems. Although the Hubbard model is simple, 
it captures main features of the electron correlations in narrow energy bands. 8,9 ) 
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In §2, we introduce the model for studying the Mott transition and see the physics of the non- 
interacting case. In §3, the previous studies on this model are reviewed. In §4, PIRG, which is a 
numerical algorithm we use in this study, is summarized. In §5, we show PIRG results and deter- 
mine the critical values of electron-electron interaction and discuss the nature of metal-insulator 
and magnetic transitions. An important conclusion from the numerical results is the existence 
of a nonmagnetic insulator phase sandwiched by a paramagnetic metal and an antiferromagnetic 
insulator. The Mott transition occurs at smaller U/t than that of the antiferromagnetic transition. 
We also discuss the physical properties near the two transitions. In particular, anomalous metallic 
state with strong renormalization of carriers is suggested near the metal-insulator transition at 
large U and large t'/t. The first-order transitions are suggested both for the antiferromagnetic 
and metal-insulator transitions at t'/t = 0.2. They are in sharp contrast with the filling-control 
transition where a concurrent single continuous transition at the vanishing doping concentration is 
observed with critical enhancement of charge and magnetic fluctuations with scaling properties. n ) 
In §6, we summarize our conclusions. 

§2. Model 

2. 1 Hamiltonian 

In this chapter, we introduce an extended Hubbard model with longer-ranged transfers for the 
study of the magnetic and metal-insulator transitions. We use the following Hamiltonian. 

H = H t + Hu 

H t = - j2 ( c h c j° + h -°- 

<i,j>,<7 
N 



tij — < 



t = 1.0 for are the nearest neighbor sites 

—t' for are the next-nearest neighbor sites (2.1) 
otherwise 

where i and j represent the lattice points, c| CT (cj CT ) the creation (annihilation) operator of an electron 
with spin a on the i-th site, rii a = cl a Ci a , tij the transfer integral between the i-th site and the j-th 
site, U the on-site Coulomb interaction and iV the number of lattice sites. We take t = 1.0 as the 
energy scale. We study two-dimensional square lattice with the next-nearest-neighbor transfer t' 
illustrated in Fig.l. 

There are many simplifications in the Hubbard model. The most drastic one is to consider only 
electrons in a single orbit. However, low-energy and low-temperature properties of a simple band 
near the Fermi level even in the presence of a complicated structure of high-energy bands far from 
the Fermi level are often well described by the Hubbard model. Therefore, the above Hubbard 
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model has some generality in discussing fundamental features of quantum phase transitions, which 
occurs at zero temperature in the strict sense. 





Fig. 1. (a)The lattice structure of the two-dimensional Hubbard model on a square lattice. (b)The lattice structure 
of the two-dimensional Hubbard model with next nearest neighbor transfers t' . 



2.2 Fermi surface 

Before we deal with the interacting Hamiltonian H = Ht + Hjj , we study the non- interacting one 
H = Hf. The energy band is given by 

Ek = —2t (cos k x a + cos k y a) + 4i' cos k x a cos k y a (2-2) 

where a is the lattice constant. Hereafter we take the length unit a = 1. In this system the Fermi 
surfaces at four choices of t' are shown in Fig. 2. These are calculated on the 200 x 200 meshed 
momentum plane. These figures show that the Fermi surfaces bend and the nesting is destroyed 
when t' increases. Note that in the range where t' is larger than about 0.6, another Fermi surface 
appears near (0, 0) as shown in the t' = 0.8 case. 

2.3 Range oft' 

Some studies 12-14 ) have been done to decide the three-band d-p model parameter by comparing a 
cluster calculation of the d-p model consisting of Cu-d x 2_ y 2 and 0-p a orbitals and the experimental 
values on the high-Tc cuprate superconductors. Some studies 15,16 ) have also been done to decide 
the parameter t' in Hamiltonian (2.1) by comparing the exact diagonalization results of the Hubbard 
model Hamiltonian (2.1) with those of the d-p model or by comparing tight-binding dispersion with 
the first-principles results. The result for if Jt is on Table I. k— (ET^-X" is an organic superconductor 
with two-dimensional anisotropy, whose lattice structure is different from Fig.l and is frequently 
mimicked by the structure illustrated in Fig. 3. In this paper we consider the t' = 0.2 and t' = 0.5 
cases in the lattice structure shown in Fig. 1(b). We refer to these examples to show that our choices 
of t' are within a realistic and sensible range as a model of several materials. However because our 
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Fig. 2. The Fermi surface of the non-interacting Hamiltonian Ht 



interest is on metal-insulator and magnetic transitions on a general ground, we do not focus on 
detailed comparison with experimental results. 



Table I. The value of next nearest transfer t' estimated for several materials 

t' material 

0.2 LSCO 

0.3 ~ 0.5 YBCO,BSCCO 

-0.8 k - (ET) 2 X [ET=BEDT-TTF,X=CuN(CN) 2 X',X'=Cl,Br,- • •] 



2-4 Quantum phase transitions 

In strongly correlated electron systems, metal-insulator transitions from the Mott insulator occur 
through three different routes. 1 ) In the extended Hubbard model described by Hamiltonian (2.1) 
they are simulated in the following: 

• Bandwidth control. 

— The bandwidth is controlled by the relative strength of electron-electron interaction U/t 
and/or the relative transfer t' jt in the Hamiltonian (2.1). 
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Fig. 3. The lattice structure of k — (ET^X, 2D-organic superconductor. 

• Band filling control. 

— In the single band Hubbard model without the uniform magnetization, the band filling 
n = may be controlled where M is the total number of particles and N is the number of 
the lattice sites. 

• Dimensionality and lattice structure control. 

— In our model the information of the dimensionality is stored in tij . Even if the dimension is 
fixed, the bond structure tu has a large influence on quantum phase transitions. 

The corresponding controls of the above three key parameters can experimentally be achieved and 
designed, for example, by doping, pressure, chemical composition and magnetic fields. Although 
the filling-control transition has been studied by the QMC quite extensively, 1 ' 4 ~ 6 ) the nature of 
the bandwidth-control transition has not been studied on the same detailed level as we discussed 
briefly in §1. In case of the perfect nesting at t' = as shown in Fig. 2, the critical value U c for the 
bandwidth-control transition is believed to be zero. It is also believed that the critical value U c is 
nonzero at t' ^ because of the destruction of the nesting. Then the bandwidth-control transition 
can only be studied at t' / 0. It is difficult, however, to apply QMC in t' ^ regime because of 
the sign problem. 

Since the bandwidth control is one of the main routes of the metal-insulator transition in many 
materials as discussed in §1, it is desired to clarify its general feature. Here we focus our study on the 
bandwidth-control transition since PIRG for the first time makes it possible to study it numerically. 
More concretely, we fix the band filling at n = 1 and study two cases of the relative transfers 
t'/t = 0.2 and 0.5 as we discussed. We search the quantum phase transition point by changing the 
relative strength of correlations U/t. We determine the critical values U c of the transitions scaled 
by t = 1.0. We study two quantum phase transitions, one the metal-insulator transition and the 
other, magnetic transition. The former is studied by estimating the U dependence of the double 
occupations and the charge excitation gap while the latter by measuring the equal-time spin-spin 
correlations. 
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§3. Previous studies on this model 

There exist some studies on the Hubbard model (2.1) in the literature. The critical values U c 
determined by these previous studies are on Table II and Table III. The Hartree-Fock approxima- 
tion is the simplest way to approach this issue. 17 ) Kondo and Moriya have used the Hartree-Fock 
approximation for the two-dimensional model and calculated the free energy as a function of the 
uniform staggered moment m = \ {n^) — {nj\)\/2 and the strength of interaction U. They have 
searched the case of U in which m 7^ gives the lowest energy, to discuss an antiferromagnetic 
transition with the periodicity given by the wavenumber (71", tt). To discuss metal-insulator transi- 
tions, they have studied the presence of the charge excitation gap in the density of states. They 
have also studied the same system by using the Green's function approach with the self-energy up 
to the second order in the strength of the interaction U. 18 ^ 

There exist also several studies on the paramagnetic state and the antiferromagnetic state by the 
Gutzwiller approximation. 19 ~ 22 ) It was claimed 17 ) that since difference in energy is small between the 
antiferromagnetic states estimated by the Hartree-Fock approximation and that by the Gutzwiller 
approximation, it would be possible to estimate the critical value of the magnetic transition by using 
the paramagnetic state of Gutzwiller approximation and the antiferromagnetic state of Hartree- 
Fock approximation. A general tendency of the Hartree-Fock approximation shows the existence of 
the antiferromagnetic metal phase. However, the Hartree-Fock and Gutzwiller approximations may 
overestimate the region of symmetry-broken states such as the antiferromagnetically ordered phase. 
Therefore, careful studies are desired to determine the phase diagram by considering quantum 
fluctuation effects. 

The method to discuss Mott transitions by QMC shares similarity to the method we use in 
PIRG study shown in the next section. Quantum Monte Carlo suffers from the sign problem in 
this regime of t'/t. As far as we know, no QMC studies for the bandwidth-control metal-insulator 
transition are available. There are some QMC studies for magnetic transitions in weakly frustrated 
regime. 5 ' 23 ^ The magnetic transition has been discussed by the equal-time spin correlations only 
for the possibility of the antiferromagnetic order at the wave number (tt, tt) while the possibility 
for other magnetic states has not been examined. The definition of equal-time spin correlation is 
shown in the next section. 

§4. Path-integral renormalization group 

Path-integral renormalization group (PIRG) is a numerical algorithm for studying the ground 
state properties. 7 ) The process projecting onto the ground state \ip g ) is performed in the imaginary 
time direction by the Hamiltonian operator H as, 

1^) = }™ e *P[- T H]\ ( t ) initiai)- 
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Table II. Critical value U c /t in case of t' /t = 0.2 





magnetic transition 


Metal-Insulator transition 


Hartree Fock 17) 


2.064 


2.064 


Green's function 18 - 1 


4.00 


4.00 


Gutzwiller 17) 


3.920 


no data 


Quantum Monte Carlo 23 ' 


2.5 ±0.25 


no data 


Table III. Critical value U c /t in case of t' /t = 0.5 




magnetic transition 


Metal-Insulator transition 


Hartree Fock 17) 


3.215 


3.290 


Gutzwiller 17) 


5.259 


no data 



Therefore its formalism can be applied to any kind of systems and there is no restriction on the 
spatial dimension of the system. In this path-integral formalism, the ground state is represented 
by a set of chosen basis states \4>i), 

i 

By the numerical renormalization, relevant basis states are selected and irrelevant basis states are 
integrated out. This makes it possible to calculate an approximate ground state within a fixed 
dimension, L, of the Hilbert space, directly as an optimized linear combination of the numerically 
chosen basis states: 

L 
i=l 

In this process, the negative sign problem does not appear even in generic fermion models. In this 
way, PIRG can be applied to the systems which can not be treated by existing algorithms such as 
quantum Monte Carlo method or the density matrix renormalization group. 

Because converged states by PIRG are an approximate ground states under the restriction on 
the number of the basis states, the properties in the exact ground state can be achieved by the 
extrapolation of the number of states L to the dimension of the whole Hilbert space. We have 
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shown a general extrapolation procedure 10 ) by using the relation: 

(H) - (H) g oc AE 

where AE is the energy variance, 

AE = ± 1 , x 1 , 

(H) 2 

( ) g , the expectation value in the true ground state and ( ), the expectation value in an approx- 
imate ground state. This relation holds for sufficiently converged approximate state. When we use 
the Slater determinants as the basis states \4>i), we confirm that the results with L roughly more 
than a hundred follow the above relation and a linear extrapolation can be used in the Hubbard 
model with the parameter values we employ in this paper. On the physical quantity A, a similar 
relation holds in most cases. 10 ) 

(A) - (A) g oc AE. 

For short-ranged correlation functions, the linearity holds at the same level as the energy in the 
thermodynamic limit. 10 ) By these extrapolation procedures, more accurate results are obtained. 
We note that, although the results at finite L satisfy the variational principle, it is not strictly 
satisfied for the extrapolated results. 

Error bars of estimated physical quantities in our method have basically two origins. One comes 
from the extrapolation procedure from finite L to the dimension of the whole Hilbert space. If 
the linearity in the fitting as a function of the energy variance becomes worse, the error bars 
increase. The other comes from the extrapolation to the thermodynamic limit by assuming finite 
size corrections. The latter is a common origin of error bars to all the other numerical methods for 
finite size systems while the former is specific to this PIRG method. The error bars in our analyses 
are given from the combination of these two types of error bars. 

§5. PIRG results 

5. 1 Method for analyzing transitions 
5.1.1 U dependence of energy 

When a transition caused by the on-site interaction U is of the first-order, dE g /dU should have 
a jump at the transition point U c because of the level crossing of the ground state as a function 
of U. Here E g is the ground state energy, namely {H) g . Even if the transition is of higher-order, 
some singularities of the energy may be seen at U c in general. From Hamiltonian (2.1), dE g /dU 
can be transformed as 



9Eg 

dU 



E n *t n *i- ( M - x)^ > ( 5 - x ) 
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where M = Mf + Mj is the particle number. Then, we can determine the transition point by 
estimating the double occupancy: 

1 N 

( n 1 n l)g = J^Y,( n ^ n n)g- ( 5 - 2 ) 

i=l 

In our calculation, the double occupancy is obtained after the extrapolation to the thermodynamic 
limit by the fitting in the form (n^ni) g (N = oo) = (n-tni) g (N = N) + a/s/N. The data we have 
obtained for the double occupancy indeed well follow this finite-size correction as in energy. 
5.1.2 Metal-insulator transitions 

Metal-insulator transitions can be discussed by computing the presence of a jump of the chemical 
potential at half filling. Because the chemical potential is defined by [i (n) = dE/dn where n = 
(Mf + Mi) /N is the particle density, we calculate the chemical potential by the following equation: 

2M-1\ _ AE _ Eg (M, M) — E g (M — 1, M — 1) 



) = " ~ < 5 - 3 > 

where E g (Mf , Mf ) is the ground state energy of the whole system including Mf electrons with the 
up spin and Mj electrons with the down spin. The charge excitation gap, if it exists at half filling, 
is calculated as 

ir / TV 4- 1 \ / A7" _ i M 

(5.4) 



N + l\ (N-l 
A* 1 » - A* 



A C = I 

2 L V N J V N 

Namely, we calculate the ground state energy E g (N/2 - 1, N/2 - 1), E g (N/2,N/2), 
Eg (N/2 + 1, N/2 + 1) by PIRG after extrapolation to the thermodynamic limit. We estimate 
the charge excitation gap from Eq.(5.4). We note here that the accuracy becomes worse for the 
estimate of the charge-excitation gap because we have to perform the second-order differentiation 
numerically from the ground state energies for the particle numbers 2(M + 1), 2M and 2(M — 1). 
We use a finite-size scaling function in the form 

A c (AT) = A c + J=A' (5.5) 

where iV is the number of sites in two-dimensional systems. We take the above scaling form 
because the finite-size effect for the Hartree-Fock SDW gap equation is given by the series of 
1/y/N. Equation (5.5) is the same as that used in the literature. 5 ^ 

In addition to the charge excitation gap, we also calculate the momentum distribution to discuss 
the electronic structure of systems. It is defined by 

TV 



™(9) = ^E (4t c *t + °W e^ R ^ (5.6) 
where Ri is the vector representing the place of the i-th site. 
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5.1.3 Magnetic transitions 

Magnetic transitions can be discussed by the equal-time spin correlations in the momentum space, 
S(q), defined by 

1 N 

s(</) = w E(^> el9M) ( 5 - 7 ) 

where S{ is the spin of the i-th site and Ri is the same as that used in Eq.(5.6). Each element of 
the spin is defined by 

sf = ^(sr+ s-) = \ (cJt^ + 4^) 

S! = \{n l ^-n ll ). (5.8) 

If a long-range magnetic order exists with the momentum q peak in the thermodynamic limit, it 
can be extracted from the extrapolation of S{q peak ) in finite-size systems. From the theory of spin 
wave in two dimensions, the following scaling is expected: 

S{q P eak) = Nm 2 + S C VN (5.9) 

where m is the staggered magnetization 

m = — ^exp(iq peafc r i )(n iT - n^), (5.10) 

i 

N, the number of lattice sites of the two-dimensional system, S c , the short-ranged part of the 
structure factor and (xi,yi), the vector indicating the coordinate of the i-th site. Then if the 
long-range order exists at q pea k, we expect S{q peak )/N to follow a linear function of l/y/N with a 
nonzero offset at 1/ y/~N = 0. 

5.2 t'/t = 0.2 case 

5.2.1 Metal-insulator transition at t'/t = 0.2 

Figure 4 shows the double occupancy as a function of the relative interaction U/t after the size 
extrapolation to the thermodynamic limit. It suggests a jump of (n-j-n^) between U/t = 3.2 and 3.3. 
Although the possibility of a very sharp second-order transition can not be excluded, a first-order 
transition seems to occur at U/t = U c ±/t = 3.25 ± 0.05. The jump of the double occupation (n-fn|) 
is estimated to be roughly between 0.04 and 0.05. 

Here, we discuss the origin of relatively large error bars near the transition points in Fig. 4. Near 
the first-order transition, fluctuations are in general not enhanced. Therefore, one might argue 
that the large error bars observed in Fig. 4 as compared to Fig. 12 would contradict this general 
expectation. In the present case, however, the large error bars are actually generated from the 
character of the first-order transition itself. The level crossing between metallic and insulating 
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phases seem to occur at different U c for different system sizes. U c seems to increase with increasing 
system sizes. This is natural because the imposed periodic boundary condition appears to suppress 
charge fluctuations in smaller system sizes. Then, if we fix U, there exists a region where the 
insulating behavior at smaller system size is transformed to a metallic state at larger system sizes. 
In fact, for example, the ground state at U = 3.0 indicating a signature of insulating behavior with 
(n^ni) ~ 0.15 at the 6 by 6 lattice transforms to a metallic behavior with (n^n^) ~ 0.2 at the 10 
by 10 lattice. This transition makes the extrapolation difficult and makes the errors large. 

Figure 5 shows the charge excitation gap near this transition. This figure suggests that the charge 
excitation gap around 0.1 exists at U/t = 3.5 and does not exist at U lower than 3.2 in agreement 
with the above analysis by dE g /dU. Therefore we conclude that the first-order transition shown 
in Fig.4 is a metal-insulator transition. A sharp drop of the gap with increasing system size for 
U/t = 3.0 and 3.2 is related with the tranformation from the insulating to the metallic branch 
mentioned above. 

Here we estimate the validity of the value of the charge excitation gap. From Eq.(5.4), the charge 
excitation gap A c is equal to the half of the sum of the excitation energies for making a holon and 
a doublon in the insulating state. Consequently, if we add a doubly occupied site in a half-filled 
system, the excitation energy is equal to 2A C . Then the total energy cost caused by the jump of 
averaged double occupancy in the system to make a metallic state from an insulating ground state 
is estimated to be roughly 2N5(n^ni)A c , where 5(n^ni) denotes the jump of (n^ni) between the 
insulating and metallic states. Since this estimate does not include the interaction energy among 
holons and doublons, this energy cost may be the lower bound for the true energy cost to make a 
metal from the insulating ground state for U > U c \. The true energy cost can be computed roughly 
as (U - U cl ) x S^ni) which is proportional to (U - U cl ) x | [dE g /dU] meta i - [dE g /dU] insulator \. 
Therefore, 

2iV5(n T n i )A c < N(U - U cl ) x 5(n^ ni } 

* A c < ^=f* (s.ii) 

From this equation, the charge excitation gap A c at U/t = 3.5 is estimated to be less than (3.5 — 
3.3)/2 = 0.1. The result shown in Fig.5 satisfies this rough estimate. It should be noted that even 
in the first-order transition, the charge excitation gap opens continuously from zero for U > U c \ as 
we see from Eq. 5.11. 

In addition to that, the behavior of the momentum distribution at half filling shown in Fig. 6 
lends another support to the above discussion. At U/t = 3.2, the momentum distribution is close 
to 1.0 inside the Fermi surface of the non-interacting system shown in Fig.2 and nearly outside. 
At U larger than 3.5 the jump smears and the ground state of the system appears to change to a 
non-Fermi-liquid state. Similar behaviors are shown for the simple square-lattice Hubbard model 
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in Fig.20 of Ref. 4 ) Because we do not know whether a metallic state which is not the Fermi liquid 
exists or not and can not determine the jump of the momentum distribution precisely, Fig.6 does 
not directly justify that the system is in an insulating phase for U/t > 3.5. However, the figure 
indicates a qualitative change in the electronic structure around U c i = 3.25 ± 0.05. 
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Fig. 4. The expectation value of the double occupancy {n^n^g after finite-size scaling. The plots are for t'/t = 0.2. 
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Fig. 5. Finite-size scaling of A c for some choices of the strength of the electron-electron interaction U in t'/t = 0.2 
lattice. 



5.2.2 Magnetic transition at t'/t = 0.2 

Examples of the structure of the equal-time spin correlations are shown in Fig.7. The peak of S(q) 
at half filling is always at (it, it) in the momentum plane. Therefore we study the antiferromagnetic 
correlation at wavelength (7r,7r) for t'/t = 0.2. 
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Fig. 6. The momentum distribution on 10 x 10 lattice at half filling for t/t' = 0.2 

The finite-size-scaling extrapolation of S(tt, tt) is shown in Fig.8. We also show finite-size scaling 
in other type of plots. S(ir, tt) goes to infinity in the thermodynamic limit when the antiferro- 
magnetic order exists. Figures 10 and 11 show these behaviors. From these finite-size scaling, the 
magnetic transition is estimated to occur at U&ft = 3.45 ± 0.05 for t'/t = 0.2. 

This value is substantially larger than TJ^/t = 2.5 ± 0.25 claimed from QMC study. 23 ) Because 
systems only up to N = 64 have been studied by QMC in that study and Fig.8 shows the sharp 
decreases of S(tt,tt) on the systems larger than 8 x 8, we think that our estimate gives better 
accuracy and the true transition point U C 2/t is at 3.45 ± 0.05 for t'/t = 0.2. A sharp decrease 
of S(tt, tt) is again related with the finite size effects accompanied with the first-order magnetic 
transition: For U/t = 3.0 and 3.2 the smaller systems show the behavior of AF ordered state, while 
it is lost in larger system sizes. The first-order character of the magnetic transition is supported 
from the following two analyses. One is the presence of a sharp decrease of the magnetization at the 
transition. Figure 9 shows the results after the finite-size scaling extrapolation for S(ir, tt)/N, which 
is equal to m 2 . This suggests either a sharp continuous transition or the behavior of the staggered 
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magnetization m at the first-order transition. The second support is in Figs. 10 and 11 which show 
that when the electron-electron interaction U increases, S(tt, it) goes to diverge suddenly between 
U/t = 3.4 and U/t = 3.5, not continuously. The data at U/t = 3.0 ,3.2 and 3.4 do not show growth 
of short-ranged correlation at the largest system size. 

Although we do not detect a jump of dE g /dU at this first-order magnetic transition in Fig.4, it 
is likely that the energy difference between paramagnetic insulator and antiferromagnetic insulator 
is far less than that between paramagnetic metal and paramagnetic insulator and then a jump of 
dEg/dU is too small to be detected in the present accuracy of calculations. 



2* U/t = 3.2 2* U/t = 3.4 




Fig. 7. The equal-time spin correlations in momentum space on 10 x 10 lattice at half filling for t'/t = 0.2 and 
U/t = 3.2, 3.4, 3.5 and 4.0. 



5.3 t'/t = 0.5 case 

5.3.1 Metal-insulator transition at t'/t = 0.5 

Figure 12 shows the double occupancy after the size extrapolation as a function of the relative 
interaction U/t. A sharp change of the double occupancy occurs between U/t = 4.5 and 5.0. 
Differently from t'/t = 0.2, this change does not seem to be a jump. Hence, this transition at 
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Fig. 8. Finite-size scaling of S (tv, n) for some choices of the strength of the electron-electron interaction U in 
half-filled t'/t = 0.2 lattice systems. 
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Fig. 9. S (tt, 7r) /N in the thermodynamic limit and as a function of the strength of the electron-electron interaction 
U in half- filled t'/t = 0.2 lattice systems. 



U c i/t = 4.75 ± 0.25 seems to be a continuous transition. 

Figure 13 shows the charge excitation gap at larger U than this transition point. Although the 
finite size scaling is difficult, we can estimate the value of the charge excitation gap as in the 
same way as we have done at t'/t = 0.2 using Eq.(5.11). According to this equation, the charge 
excitation gap A c at U/t = 6.0 should be comparable or less than (6.0 — 5.0)/2 = 0.5 and A c shown 
in Fig. 13 satisfies this prediction. Hence judging from the rough U dependence of the charge gap, 
we conclude that the transition seen in Fig. 12 is the metal-insulator transition. 

We also show the momentum distribution in Fig. 14. Although the qualitative change of the 
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Fig. 10. Finite-size scaling of 5 (tt, it) for some choices of the strength of the electron-electron interaction U in 
half-filled t'/t = 0.2 lattice systems. 
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Fig. 11. Finite-size scaling of S (ir, ir) for some choices of the strength of the electron-electron interaction U in 
half-filled t'/t = 0.2 lattice systems. 



momentum distribution at U/t = U&/t = 4.75±0.25 can be seen, it is remarkable that the change at 
the metal- insulator transition point is smaller and slower as a function of U/t than the case at t'/t = 
0.2. It may be related to the order of the transition; the fluctuation at the continuous transition is 
larger than that at the first-order transition and makes the change gradual. The behavior of the 
momentum distribution at U/t = 4.5, t' /t = 0.5 implies that a strong renormalization of carriers 
may grow already in the metallic state near U c \ and an anomalous metal state may be stabilized. 
Spin fluctuations are not critically enhanced in this region because the spin correlations are short- 
ranged even in the insulating side near the phase boundary. Therefore, unusual metallic behavior 
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is not ascribed to the spin fluctuations. The origin of the strong renormalization may come from 
solely charge fluctuations near the metal-insulator transition. Superconducting pairing correlations 
are intriguing quantities to be studied in the future in terms of the instability of this unusual metal. 
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Fig. 12. The expectation value of the double occupancy {n-f n{)g after finite-size scaling. The plots are for t' /t = 0.5. 
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Fig. 13. Finite-size scaling of A c for some choices of the strength of the electron-electron interaction U in t' /t — 0.5 
lattice. 



5.3.2 Magnetic transition at t' jt = 0.5 

The structures of the equal-time spin correlations S{q) are different from those for t'/t = 0.2. 
Although the systems are at half filling, Fig. 15 shows incommensurate peaks near (it, it) for U/t < 
7.0. The equal-time spin correlations have a single peak at (ir, it) for 6 x 6 and 8x8 lattices and 
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Fig. 14. The momentum distribution on 10 x 10 lattice at half filling for t' /t = 0.5 



have incommensurate peaks for 10 x 10 and 12 x 12 lattices. We confirm that this tendency holds 
for U/t = 4.5, 5.0, 5.5, 6.0, 6.5 and 7.0. 

First we analyze the incommensurate peaks. The incommensurate peaks shown in Fig. 15 have 
asymmetric momentum dependence in terms of the four fold symmetry in x or y direction. These 
do not, however, have a physical significance of the symmetry breaking because the system size 
is finite. The origin of this asymmetry may be that PIRG can treat only a subspace of the total 
Hilbert space and even after the extrapolation L — ► oo, this tendency remains. We discuss this 
problem below. 

The similar incommensurate peaks appear away from half filling for the square lattice with 
t'/t = 0. 5 ) In that case, the perfect nesting is destroyed by the hole doping and some other 
partial nesting occurs instead. In the present case, however, the perfect nesting is destroyed by the 
frustrating term t' and no clear nesting vectors at the Fermi level are found in the non-interacting 
system. 

To discuss the presence of the long-range order, we have to study these superimposed incommen- 
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surate peaks in the thermodynamic limit. In the following analysis, we assume that the occurrence 
of the incommensurate peaks is related to the flat structure of the peak at (tt, tt) of the equal-time 
spin correlations for the non-interacting system as shown in Fig. 16. Though for 8x8, S(tt, it) shows 
a single highest peak, for 10 x 10 and 12 x 12, momenta q at which the S(q) has the highest value 
are (it, it) as well as the nearest momenta to it. This is a similar system-size dependence to that 
of the incommensurate peaks for the interacting system. Although for 4 x 4, 6 x 6 and 8x8, there 
are single peaks at (tt,tt), for 10 x 10 and 12 x 12, the incommensurate peaks as shown in Fig. 15 
appear at the nearest momenta to (tt,tt). We also confirm that the 'plateau' behavior holds for 
30 x 30 on non-interacting system. Therefore the incommensurate peaks are likely to appear also 
at the nearest momenta to (tt, it) for the interacting systems larger than 12 x 12 and the momenta 
at which peaks appear can be speculated to converge to (it, tt) in the thermodynamic limit. For 
this reason we just illustrate the peak values for S(q) on the same plane in Fig. 17. From this figure, 
it is clear that the behaviors of the equal-time spin correlations are different between U > 7.5 
and U < 7.0. Although the commensurate peaks of the equal-time spin correlations grow when 
the system size increases for U > 7.5, those for U < 7.0 remains finite clearly. We conclude that 
the incommensurate peaks which appear for U < 7.0 on the 10 x 10 and 12 x 12 lattices remain 
short-ranged in the thermodynamic limit and may converge to the commensurate peak. 

As we discuss in § 5.5, the asymmetry of the incommensurate peaks may come from a precursor of 
the four-sublittice order with Bragg peaks at (tt, 0) and (0, tt). When the short-ranged fluctuations 
around (tt, 0) and (0, tt) occur in the Hilbert space far apart each other, PIRG may reproduce only 
one part of this fluctuation because of a limited number of L. Since PIRG may be regarded as 
a systematic improvement of the Hartree-Fock approximation, and the Hartree-Fock calculation 
always overestimate the ordered phase, the absence of the long-ranged order observed by PIRG 
should be reliabe for this incommensurate fluctuations. 

Next we discuss the commensurate long-range order. The incommensurate peaks disappear and 
the commensurate antiferromagnetic peak at (it, it) appears for U > 7.5 as we see in Fig. 15. To 
discuss the long-range order, we use the finite-size scaling function Eq.(5.9) again and Fig.18 shows 
the result. This figure shows that the commensurate peak which appears for U > 7.5 indicates 
the long-ranged order. The transition bears continuous character. The reason for that we can not 
detect a jump of dE g /dU is supposed to be the same as the reason for the case t'/t = 0.2; the jump 
is too small. 

5.4 Comparison with Hartree-Fock results 17 ^ 

The critical values of the relative interaction U c /t determined above is shown in Tables IV and V. 
The magnetic transition appears to occur at larger U than that for the metal-insulator transition. 
The paramagnetic insulating state seems to appear in between. This is in contrast with the Hartree- 
Fock calculation where the two transitions occur concurrently at t'/t = 0.2 while the two transitions 
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2* U/t = 6.5 2 * U/t = 7.0 




Fig. 15. The equal-time spin correlations in the momentum space on 12 x 12 and t' /t = 0.5 lattice at half filling for 
U = 6.5, 7.0, 7.5 and 9.0. Note that the asymmetry is caused by the PIRG. 



Table IV. Critical value U c /t for magnetic transitions to a commensurate antiferromagnetic state and metal- 
insulator transitions at half filling for t' /t = 0.2. 





magnetic U C 2/t 


metal-insulator U c i/t 


PIRG 


3.45 ± 0.05 


3.25 ±0.05 


Hartree Fock 17 ' 


2.064 


2.064 


Green's function 18 ' 


4.00 


4.00 


Gutzwiller 17) 


3.920 


no data 


Quantum Monte Carlo 23 ' 


2.5 ±0.25 


no data 
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U/t = 0.0 
t'/t = 0.5 




U/t = 0.0 
t'/t = 0.5 




U/t = 0.0 
t'/t = 0.5 




U/t = 0.0 
t'/t = 0.5 




Fig. 16. The equal-time spin correlations in the momentum space for t'/t = 0.5 at half filling for the non-interacting 
system. The systems are 8 x 8, 10 x 10, 12 x 12 and 30 x 30 lattices. For lattices equal to or larger than 10 x 10, 
the peak always has five-fold degeneracy around (7T,7r) which leads to the 'plateau' structure around (tt,tt). 
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Fig. 17. Finite-size scaling of S (q peak ) on some cases of the strength of the electron-electron interaction U in 
t'/t = 0.5 lattice. 
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Fig. 19. S (ir, it) /N in the thermodynamic limit as a function of the strength of the electron-electron interaction U 
in half- filled systems at t'/t = 0.5. 



occur in the opposite order and an antiferromagnetic metal is stabilized at t'/t = 0.5. This may 
be interpreted from the fact that the magnetic long-ranged order is more strongly overestimated 
in the Hartree-Fock approximation than the metal-insulator transition because fluctuation effects 
are more serious for the antiferromagnetic symmetry breaking. 

The Hartree-Fock approximation also shows that the combined transition at t'/t = 0.2 is of 
the first order. 17 ) PIRG results appear to show that the two transitions are both of first order at 
t'/t = 0.2. At t'/t = 0.5, the Hartree-Fock approximation study claims that the magnetic transition 
is of second order and the PIRG results does not seem to contradict this. PIRG also shows that 
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Table V. Critical value U c /t for magnetic transitions to a commensurate antiferromagnetic state and metal- insulator 
transitions at half filling for t'/t = 0.5. 



magnetic metal-insulator 



PIRG 7.25 ± 0.25 4.75 ± 0.25 

Hartree Fock 17) 3.2f5 3.290 

Gutzwiller 17) 5.259 no data 



the metal-insulator transition seems to become continuous. 

5.5 Presence of spin liquid insulator near the metal-insulator phase boundary 

According to the above discussion, there is a nonmagnetic insulator phase between metal-insulator 
and magnetic transitions. We discuss this remarkable fact in detail for the case of t'/t = 0.5. There 
have been not many studies on the Hubbard model on a square lattice with nearest-neighbor and 
next-nearest-neighbor transfers. However, extensive studies have been done on spin-1/2 Heisenberg 
model on the same lattice 24-29 ) which is called J1-J2 model, although the result is not sufficiently 
conclusive due to the negative sign problem and lack of efficient numerical algorithm. The Hamil- 
tonian of the J1-J2 model is 

H = "A? ^ • S j 

ij 

J\ > for i and j are the nearest neighbor sites 

J2 > for i and j are the next nearest neighbor sites (5-12) 
otherwise 



Jij 



where Sj is the quantum spin operator (S = rj) of the i-th lattice site. The ground state of the 
classical J1-J2 model on the square lattice is the two-sublattice Neel state for J2/J1 < 0.5 and 
the degenerate four-sublattice Neel state for J2/J1 > 0.5. Similarly to the Hubbard models, the 
quantum model with J\ only has the Neel antiferromagnetic long-range order in the ground state. 
According to the literature, 28 ' 29 ) increasing the frustration J2 destabilizes the Neel order and at 
J2/J1 — 0.4 a phase transition is inferred to occur to some different state and at J2/J1 — 0.6 a 
phase transition to a collinear order, which is a special case of the four-sublattice Neel state, occurs. 
Although no conclusive results are available, the state stabilized for 0.4 < J2/J1 < 0.6 is speculated 
to be a nonmagnetic (or spin liquid) insulator. In fact, a translational symmetry breaking due to 
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a dimer order has been suggested, 30 ' 31 ) although it is controversial. 32 ) 



It is known that the half-filled Hubbard model can be treated by the second-order perturbation 
and in the limit of large interactions, the Hubbard model is mapped to the Heisenberg model with 
the superexchange coupling, 



Now we focus on the case t'/t = OM. This parameter value leads to the mapping: 



In this region, the ground state of the J\- J2 model seems to have an antiferromagnetic long-range 
order of the Neel state. This is also supported by our result for U/t > 7.5. According to our 
result on the Hubbard model, however, there is an insulating phase without the antiferromagnetic 
long-range order for 5.0 < U/t < 7.5. The reason for the appearance of this intermediate state for 
5.0 < U/t < 7.5 in the Hubbard model is speculated in the following: In the Heisenberg model, the 
driving force to destabilize the commensurate antiferromagnetic long-range order is the interplay 
of quantum mechanical fluctuations and the frustration J2. However, in addition to this interplay, 
charge fluctuations ascribed to a finite amplitude of double occupancy reduces the magnetic long- 
range order in the Hubbard model at finite U and enhances quantum fluctuations. Then it is likely 
that the nonmagnetic insulator may be stabilized even for t' smaller than 0.63 (44> J2/J1 — 0.4) if 
U decreases from infinity. Moreover, incommensurate peaks and their symmetry breakings to x or 
to y direction may be understood in the following: S(q) of two-sublattice Neel state stabilized at 
J2/J1 < 0.4 has a Bragg peak at q = (tt,tt) and S(q) of four-sublattice collinear state stabilized 
at J2/J1 > 0.6 has a Bragg peak at q = (tt, 0) or (0, 7r). The dimer order or its fluctuations 
may also induce a peak around q = (tt, 0) or (0,7r). Hence, S(q) of a nonmagnetic insulator 
state, stabilized in the intermediate region, can be speculated to have peaks between (tt, tt) and 
(tt, 0), or between (tt,tt) and (0,7r). Incommensurate peaks shown in Fig. 15 is indeed found in this 
momentum region. Though these incommensurate peaks may appear only for finite-size systems, 
the physics contained in it may be related to a precursor and fluctuations of the collinear or the 
dimer state. Even when the dimer or plaquet orderings leading to the translational symmetry 
breakings is stabilized in the Heisenberg limit U — > 00, such translational symmetry breakings 
would be destroyed under sufficient charge fluctuations realized with reducing U near the metal- 
insulator transition. The nature of the nonmagnetic insulator phase found between the magnetic 
and metal-insulator transitions is an outstanding and fundamental problem to be clarified. This 
interesting issue will be studied in separate work. 

5.6 Phase diagram 

By using the PIRG results, the result of the frustrated Heisenberg model, and the fact that 
the system is in a paramagnetic- metal phase at U/t = while in an antiferromagnetic-insulator 





(5.14) 
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phase at t'/t = 0, we show a schematic phase diagram in Fig.20. For smaller t'/t up to 0.2, 
the metal-insulator and antiferromagnetic transitions occur more or less concurrently with a first- 
order nature similarly to the Hartree-Fock results. On the other hand, with increasing t'/t, the 
phase diagram shows qualitative difference from the Hartree-Fock prediction. The appearance of 
the nonmagnetic insulator rather has a consistency with the strong coupling Heisenberg picture. 
The present numerical results clarify how these two contradicting pictures are reconciled at finite 
U. The emergence of unusual metallic behavior observed near the phase boundary between the 
nonmagnetic insulator and the paramagnetic metal is one of the most remarkable findings in this 
phase diagram. We have not studied the spin excitation gap. This is left for future studies. 

In the phase diagram, we discuss the regimes of the organic materials and the cuprate super- 
conductors. Although the organic materials such as BEDT-TTF compounds have different lattice 
structure from the present t-t'-U model, these compounds seem to be located in the regime of 
smaller U/t and t'/t near the concurrent transition boundary, because they show concurrent tran- 
tions in many cases, while the nonmagnetic insulator except the spin-Peierls phase has not been 
observed so far. The cuprate superconductors may be mimicked by the present model while they 
undergo the filling-control transition contrary to the present bandwidth-control study. In any case, 
these superconductors with higher transition temperatures such as Bi, Tl, and Hg compounds seem 
to have the effective t'/t ~ 0.3 - 0.4 with U/t ~ 8. 16 ) This is deeply inside the AFI phase at half 
filling, namely in the mother materials. 




Fig. 20. A schematic phase diagram in the plane of U/t and t'/t. The symbols x represent points where we have 
analyzed in detail by the PIRG calculation. Antiferromagnetic insulator, spin liquid insulator (or nonmagnetic 
insulator), and paramagnetic metal are denoted by AFI, SLI, and PM, respectively. The arrow indicates the phase 
boundary between AFI and SLI in the limit U/t — > oo roughly speculated from the frustrated Heisenberg model. 
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§6. Summary and Discussion 

Mott transitions in strongly correlated two-dimensional systems may occur through various dif- 
ferent routes such as bandwidth control, band filling control and dimensionality control. The filling- 
control transition has already been studied quite extensively in two dimensions. 1 ' 4 ' 5 ) Because the 
Mott transition in the half-filled square lattice with the nearest-neighbor transfer is expected to 
occur at U = owing to the perfect nesting, it is necessary to treat the systems with longer-ranged 
transfer for studying the bandwidth-control transition. Before PIRG, no powerful numerical meth- 
ods had been available except in one-dimensional systems for discussing such a system because of 
the sign problem in the Monte Carlo method and for this reason the bandwidth-control transition 
has not been studied on the same level as the filling control case. In this paper, by applying PIRG, 
we have for the first time studied the bandwidth-control transition systematically. We have stud- 
ied the two-dimensional Hubbard model on the square lattice with next-nearest-neighbor transfers 
if to clarify the nature of the bandwidth-control transition. More concretely, we determine the 
critical strengths of the on-site Coulomb interaction U c /t of the magnetic transition and of the 
metal-insulator transition for two choices of the parameters t'/t = 0.2 and 0.5 for the Hamiltonian 
(2.1). For these systems, there exist some results from the Hartree-Fock approximation in the 
literature, although they may be uncontrolled owing to the neglect of fluctuations. As easily ex- 
pected, because of the fluctuations, the critical strengths of the interactions estimated by PIRG are 
larger than those estimated by the Hartree-Fock approximation on metal-insulator and magnetic 
transitions. The present PIRG results show that the metal-insulator transition occurs at U c \ which 
is lower than the antiferromagnetic transition point U C 2, in sharp contrast with the Hartree-Fock 
prediction. The Hartree-Fock result gives the antiferromagnetic transition at lower U than the 
metal-insulator transition and an antiferromagnetic metal appears in between the two transitions. 
On the contrary, in the region U c i < U < U C 2, the nonmagnetic insulating state appears to be 
stabilized in our results. 

For the case t'/t = 0.2, the critical value of the metal-insulator transition is U c \/t = 3.25±0.05 and 
that of the magnetic transition is U C 2 = 3.45 ± 0.05. It appears to show two separate transitions in 
our accuracy. We have analyzed the order of the magnetic transition for t'/t = 0.2. S(ir, ir) appears 
to diverge suddenly between U/t = 3.4 and U/t = 3.5 without appreciable growth of short-ranged 
correlations for U < 3.4. The transition is consistent with that of the first order as a function of 
the electron-electron interaction U as the same as the result of the Hartree-Fock approximation. 
The staggered magnetization m in the thermodynamic limit seems to have a sharp increase at the 
critical strength of the electron-electron interaction in Fig. 9 and it does not contradict the above 
conclusion, although a sharp continuous transition is not excluded. We have also analyzed the 
order of metal-insulator transition. It also appears to show a first-order transition with a jump in 
the average of the double occupation. 
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For the case t'/t = 0.5, the metal-insulator transition occurs at U c i/t = 4.75 ±0.25 and the mag- 
netic transition does at U&ft = 7.25±0.25. The difference between two transition points, \U c i— 
become wider with increasing t'. Incommensurate peaks of the equal-time spin correlation appear 
near (tv,tt) in finite-size systems for U/t < 7.0. We have analyzed these incommensurate peaks by 
analyzing the amplitude of the highest peaks. This analysis indicates that these incommensurate 
peaks do not have a long-range order. S(ir,ir) between U/t = 7.0 and U/t = 7.5 suggests that the 
magnetic transition is of more continuous type. The staggered magnetization also supports this. 
Figure 12 suggests that the average of double occupation change continuously at the transition 
point and metal-insulator transition also behaves rather as a continuous one. The behavior of the 
momentum distribution supports this conclusion by indicating the presence of strong fluctuations 
and renormalizations in the metallic phase. It suggests that some anomalous metal is stabilized at 
U/t less than 4.75. This unusual metallic behavior is one of the most remarkable feature in our 
phase diagram. 

For both cases of t', a nonmagnetic insulator phase is stabilized between the metal-insulator and 
magnetic transitions. The antiferromagnetic long-range order is stabilized for the strong interaction 
regime in accordance with the expectation from the corresponding J1-J2 Heisenberg model derived 
from the strong coupling expansion in terms of t/U and t'/U, while a spin-liquid insulator may be 
realized in this intermediate interaction regime. According to the second-order perturbation in the 
limit of large interactions, the half-filled Hubbard model with the interaction U and transfer t L j 
corresponds to the Heisenberg model with the superexchange coupling = A\tij\ 2 /U. Therefore 
the Hubbard model with t'/t = 0.5 in the strong interaction limit corresponds to the J1-J2 model 
with J2/J1 — 0.25, the ground state of which seems to have an antiferromagnetic long-range order of 
the Neel state. Although we have seen that the commensurate antiferromagnetic long-range order 
appears in the Hubbard model with t'/t = 0.5, U/t > 7.25 ± 0.25, the long-range order becomes 
absent for the weaker interaction. We interpret it in the following way: The double occupancy 
allowed at finite U enhances the quantum fluctuation and reduce the magnetic local moment. This 
charge fluctuation reduces the magnetic long-range order and antiferromagnetic order is destroyed 
before the transition to a metal. The overall phase diagram also implies that the nonmagnetic 
insulator phase in the J1-J2 Heisenberg model is continuously connected with the nonmagnetic 
insulator phase in the present results, although we reserve the possibility that the dimer ordering 
speculated in the Heisenberg limit may be destroyed again by the charge fluctuation. 

Several discussions and comments on open issues may now be in order. We note the possibility 
of observing the nonmagnetic insulating phase in real materials. As we see in the phase diagram, 
Fig. 20, the width of this phase as a function of U/t becomes very narrow for small U c and could 
easily become a concurrent transition if the coupling to lattice degrees of freedom would be strong 
or effects of orbital degeneracy is present. Several organic materials may belong to this class, 
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although further careful analysis is desired. It is also desired to explore possibility of materials 
which have larger effective U. An interesting question to be examined experimentally is whether 
the nonmagnetic insulator is generically stabilized near the metal-insulator transition point when 
the 'frustration' effect such as t' is large. The nature of the unusual metal near the nonmagnetic 
insulator is an intriguing open question to be answered in future studies. This is particulary 
interesting because the spin fluctuation is not the origin of the strong renormalization. We have 
not studied the spin excitation gap. The presence or the absence of the spin excitation gap in the 
phase diagram is left for further studies. The nature of the nonmagnetic insulator state should be 
clarified in more detail in future. 

The bandwidth-control transitions in two dimensions clarified in this paper show sharp contrasts 
with the filling-control transitions. The first clear difference is seen in the separation of the metal- 
insulator and antiferromagnetic transitions. In the filling-control case, both transitions seem to 
occur simultaneously at the doping concentration 5 = 0, if the antiferromagnetic order exists in the 
insulating phase. 5 ) However, the two transitions occur at different points when the bandwidth is 
controlled. In the filling-control transition, the antiferromagnetic short-ranged correlation critically 
grows in the paramagnetic metal with decreasing 6 as S(tt, it) oc S^ 1 . Such critical enhancement 
with continuous character of the transition is not clearly observed in the present bandwidth control 
and the first-order transition is plausible for small t'. The metal-insulator transition is also clearly 
different each other: The filling-control transition shows unusual but continuous decrease of the 
Drude weight 33 ) indicating the continuous transition with specific scaling laws satisfied, 11 ) while 
the bandwidth control shows the level crossing with the first-order transition. The proximity of the 
Mott insulator in the metallic phase with an anomalous metallic state is an outstanding feature 
of the filling-control transition while such fluctuations may not be expected for the first-order 
transition inferred in the bandwidth control transition at smaller U c /t. However, it should also 
be noted that anomalous metallic states are suggested near the transition even in the bandwidth 
control as discussed above, if the transition occurs at large U owing to large t'/t with a continuous 
character of the transition as also suggested in the momentum distribution at t'/t = 0.5. 

It would be an intriguing future problem to examine the generality of the present remarkable 
features of the bandwidth-control matal-insulator transitions in other lattice structures or in three- 
dimensional systems. The effect of orbital degeneracy is also an interesting subject to be studied. In 
real materials, there exist few examples of antiferromagnetic metals such as NiS2-xSe x and V2_ y 03, 
where the order of the transitions between the metal-insulator and antiferromagnetic transitions 
is opposite to the present case. From the present study, the existence of the antiferromagnetic 
metal seems to be ascribed to a subtlety of the higher dimensionality and/or presence of orbital 
degeneracy. 
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